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ABSTRACT 

The  interaction  of  a directed  beam  of  radiation  with  an 
absorber  is  described  using  a two  level  model;  transmitted  light 
intensities  and  level  populations  are  treated  as  functions  of 
temporal  and  spatial  independent  variables;  and  the  general 
rate  equations  are  given  as  partial  differential  equations. 
Additionally,  combinations  of  material  constants  are  identified 
as  scale  parameters  and  used  to  express  the  rate  equations  in 
nondimens ional  form  as  a hyperbolic  partial  differential  equation 
in  cononical  form.  One  integration  yields  a time* independent 
relationship  between  two  transmitted  beam  intensities  at  different 
wavelengths.  Also  described  are  several  calibration  experiments 
to  evaluate  material  parameters. 
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INTRODUCTION 


This  memorandum  describes  our  work  in  investigating  the 
performance  of  an  optically  excited  absorption  system.  A general 
two  level  model  is  vised  to  describe  the  process  of  photo- induced 
transitions  in  an  absorbing  material.  One  particular  application 
of  interest  concerns  the  use  of  photochromic  and  liquid  crystal- 
photoconductor  devices  as  unconventional  photographic  media  for 
optical  data  processing  applications  --  e.g.,  input  data  recording 
and  adaptive  spatial  matched  filters.  Since  these  devices  do  not 
require  image  development  after  exposure,  they  are  attractive 
alternates  to  film  recording  in  those  applications  where  speed 
of  recording  is  a factor.  In  addition,  a variety  of  photochomics 
find  use  in  opthalmic  components,  graphic  arts,  and  displays 
(Ref.  1) ; thus  the  results  of  this  investigation  are  also  of  value 
in  these  areas. 

In  general  the  rate  equations  descriptive  of  the  two  level 
model  have  a broad  scope  and  are  useful  as  well  in  such  other 
areas  as  laser  beam  propagation,  optical  cross-section  measurements 
(Ref.  2),  photochemistry  (Ref.  3),  atomic  and  molecular  beam 
spectrocopy  (Ref.  4) , dye  laser  excitation  (Ref.  5) , and  laser 
isotope  enrichment  (Ref.  6).  In  the  first  two  cases  interest  is 
primarily  focused  on  characteristics  of  the  attenuated  light  beam, 
for  the  latter  examples  energy  level  population  dynamics  are  of 
concern. 

This  memorandum  focuses  on  a two  level  model  representing  a 
photo-absorption  process  and  considers  spatial  as  well  as  temporal 
dependencies.  The  inclusion  of  an  independent  spatial  variable  is 
necessary  to  account  for  the  fact  that  all  absorbing  atoms  or 
molecules  do  not  experience  the  same  light  level.  Light  has 


already  been  absorbed  In  passage  through  the  material.  Furthermore, 
since  the  absorption  interaction  is  a nonlinear  process,  one 
cannot  simply  scale  over  the  different  spatial  or  time  intervals 
nor  scale  from  low  to  high  power  excitation  levels  as  would  be 
the  case  for  a linear  system.  Thus,  both  the  energy  level 
populations  and  transmitted  beam  irradiance  are  functions  of  sample 
thickness  as  well  as  time,  and  the  material  will  not  exhibit  simple 
exponential  absorption  properties.  For  example,  the  transmitted 
beam  will  initially  be  strongly  attenuated  due  to  a relatively 
large  ground  state  population  (for  a ground  state  transition). 
Subsequently  the  transmitted  beam  irradiance  could  then  increase 
due  to  changes  in  the  energy  level  populations,  simulating  a self- 
induced  transparency  effect. 

The  interaction  with  two  light  beams  at  different  wavelengths 
is  also  a topic  of  interest  since  the  spectral  absorpiton  of  the 
material  influencing  the  response  at  one  wavelength  can  affect 
the  absorption  properties  at  the  other  wavelength.  This  feature 
provides  a basic  reason  for  using  such  material  in  the  above 
mentioned  optical  processing  applications.  As  a result,  data 
recording  can  be  accomplished  within  an  insignificant  development 
time.  With  such  purpose  in  mind,  the  following  model  analysis 
could  lend  itself  towards  optimization  in  the  choice  of  materials 
and  also  suggests  a favorable  means  to  tailor  or  synthesize  new 
material  parameters  for  these  applications  (Ref.  7). 

In  brief,  this  work  deals  with  a set  of  fundamental  rate 
equations  that  encompasses  several  special  cases  previously  studied 
by  other  investigators.  A number  of  material  constant  groupings 
are  found  in  the  restructure  of  the  rate  equations  in  terms  of 
nondlmensional  variables.  In  effect,  a smaller  number  of  basic 
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(material  constants)  factors  are  used  to  describe  the  behavior 
of  the  absorber* radiation  interaction  and  to  greatly  simplify  the 
analysis.  When  two  light  beams  of  different  wavelengths  Interact 
with  the  absorber,  a time- independent  relationship  has  been 
found  between  the  beam  intensities  after  passage  through  the 
material.  This  property  suggests  that  measurements  of  light 
intensity  at  one  wavelength  can  be  used  as  a probe  to  infer  the 
light  intensity  level  at  the  other  wavelength.  In  addition,  several 
methods  are  demonstrated  that  provide  estimates  for  a number  of 
these  material  constants  from  calibration  experiments.  A photo- 
chromic  material  is  used  as  a convenient  "two-level"  absorber 
to  illustrate  the  calibration  procedures.  No  attempt  is  mad<> 
here  to  integrate  numerically  the  restructed  general  rate  equations 
for  particular  situations. 


TWO  LEVEL  RATE  EQUATIONS 

TWo  unidirectional  beams  of  excitation  and  of  probing  radiation 
are  at  normal  incidence  to  the  boundary  face  of  the  absorber 
material,  as  shown  in  Figure  1.  It  is  assumed  that  the  probe 
beam  is  much  weaker  than  the  excitation  beam.  For  generality,  the 
wavelengths  of  the  excitation  and  probing  beams,  and 
respectively,  are  considered  to  be  different.  (Both  beams  need 
not  be  of  different  wavelengths;  for  example,  the  polarization 
state  of  the  weaker  probe  beam  could  be  different  from  the 
polarization  state  of  the  excitation  beam,  assuming  that  the  absorber 
material  is  polarization  insensitive.)  It  is  also  assumed  that 
light  scattered  out  of  either  beam  is  not  significant  for  homogeneous 
materials  and  can  be  neglected.  Thus,  only  the  absorption 
properties  in  the  direction  of  propagation  are  of  concern.  For 
the  one-dimens ional  geometry  considered  here,  both  beams  are  of 
uniform  intensity  in  any  plane  normal  to  the  direction  of 
propagation  along  the  x-axis,  but  vary  with  distance  x and  time  t. 

Let  1^  » 1^  (x,t)  represent  the  excitation  light  intensity 
and  I2  ■ I2  (x,t)  the  probe  light  intensity.  If  N^  (x,t)  denotes 
the  number  of  absorber  molecules  in  the  lower  energy  level  per  unit 
volume,  and  Ny  (x,t)  the  number  density  of  absorber  molecules 
in  the  upper  energy  level,  then  the  total  number  density  is 

Nq  - NL  + Ny  » const.  (1) 

The  rate  of  decrease  of  excitation  light  intensity  in  the 
direction  of  propagation  is  given  by 

dl, 

5T  - • (t,Ll  NL  + °01  + a>  *1  <2> 
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where  ctl1,  are  the  absorption  cross-sections  at  wavelength 

for  the  lower  and  upper  levels,  respectively,  and  a represents  the 

* 

absorption  rate  due  to  a host  material  interaction,  if  any. 
Similarly  the  rate  of  decrease  of  probe  light  intensity  is  given 
by 


dl 


5T  - - <”l2  \ + aU2  «B  + h 


(3) 


where  the  absorption  constants  have  similar  meaning  to  those  for 
Eq.  (2). 

Changes  in  the  number  densities  NL  and  are  given  by  the 
population  rate  equations  as 


St  * 0L1  *L1  NL  I1  ' aUl  *U1  NU  I1 


+ aL2  *L2  NL  l2  ” ctU2  *U2  NU  I2  " NU/t  (4a) 


and 


(4b) 


The  quantities  *Ll,  4>u;l,  <>u2  are  the  quantum  efficiencies  for 

each  particular  absorption  reaction  and  t is  the  "life  time"  or 
decay  time  constant  for  the  upper  energy  state.  As  a realistic 
and  practical  situation  we  shall  assume  that  both  excitation  and 
probe  beams  are  uniform  across  the  plane  of  incidence  and  constant 
in  time.  Consequently,  the  initial  and  boundary  conditions  are 

* 2 

The  units  for  1^  and  I2  are  photons /secj cm  . The  factor  hv  is 

used  to  convert  to  units  of  ergs/sec-cm  or  watts /cm  , when  h 

is  Planck's  constant  and  v is  the  frequency  of  light. 
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Ix  <0,t) 

I2  (0,t) 

(k,0) 

Nl  (x,0) 

In  the  event  that  the  material  is  initially  unexcited  then  NyQ  — 0 
auo  - No. 

Other  investigators  have  worked  with  less  general  two  level 
models.  Variations  £rom  their  results  are  due  to  specific  assumptions 
made  concerning  the  importance  of  selected  terms  in  the  governing 
equations.  It  is  worth  while  to  note  briefly  these  differences 
from  the  more  general  rate  equations  stated  above.  In  dealing  with 
a photochromic  material  Baldwin  (Ref. 8)  treats  changes  in  the 
upper  energy  level  number  density  due  only  to  absorption  transitions 
from  the  lower  to  the  upper  level  and  spontaneous  decay  down  from  the 
upper  level.  This  is  equivalent  to  setting  %i  * °l2  *l2  * 

°U2  <*>U2  * 0 in  EtJ*  (4).  In  addition,  he  also  takes  a ■ 0 in  Eq.  (2) 
and  aL2  « P ■ 0 in  Eq.  (3).  The  photochromic  rate  equations  studied 
by  Mohn  (Ref.  9)  are  equivalent  to  setting  * 0 in  Eq.  (4) 

and  gL2  ■ P ■ 0 in  Eq.  (3).  Tomlinson  (Ref.  10)  deals  with  a 
model  set  of  equations  equivalent  to  setting  1/t  - 0 and  does  not 
include  the  influence  of  a probe  irradiation;  i.e.,  terms  containing 
12  and  Eq.  (3)  are  not  included  in  his  analysis.  Together  with  a 
clever  transformation  of  variables,  he  demonstrates  an  integral 
relationship  designated  as  the  "Photochromic  Function."  In 
Ohtsuka' s (Ref.  11)  formulation  of  the  rate  equations,  pertaining 
to  laser  beam  propagation,  there  is  no  distinction  made  between 
absorption  coefficients  since  both  probe  and  excitation  wavelengths 
are  the  same.  In  addition,  his  analysis  assumes  the  condition 
equivalent  to  setting  cr^  0^. 
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It  is  advantageous  to  introduce  other  variables  and  parameters 
in  the  rate  equations.  Let 

Jx  - IL  eax,  Jx  (0,t)  - Ax,  J2  - I2  ePx,  J2  (0,t)  » A2, 


q - Ny/^,  T]  (y,0)  - 0,y  - x 


Equations  (2),  (3),  and  (4a)  become 


' ,1+  ft  -1)  'h 


fu  + .gU2l°U 

aLl  aLl 


) 


Jo  (8) 


dl  " aLl  ^Ll  <1"T»>  J1  6XP  ("ay/NoaLl)  ' 0U1  ^Ul  ^ J1  eXp  (“ay/NoaLl) 


+ aL2  ^L2  (1"T))  J2  6XP  ('Py/NoCTLl)  ’ °U2  *U2  n J2  exp  (“Py/NoaLl) 


- n/T 

With  the  further  substitutions 


J10  * 1//t  ctL1  ^Ll  ’ J20  “ 1^T  ctL2  *L2  * * " 


" J^/ JlO»  " A, /Jin,  s » s (0, A)  ™ A^/J 


TJ10 


2 20' 


2 20 


aUl  *U1^CTL1  *L1’  *2  " aU2  *U2^aL2  *L2’ 

^Ul^Ll*  1’  a2  " aL2/aLl’  S3  " ^°U2  ” aL2^aLl 
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Equations  (7)  through  (9)  are  reduced  to 


t£  - * <l  + al  1)  p 

(ID 

- • <«2  + a3  1>  8 

(12) 

>))  p exp  (-cy/N^j) 

+ (l-q-Kj  q)  s exp  (^y/H0<^j)  - q 

(13) 

Further  simplifications  can  be  made.  Introduce 
variable  q(y,X)  from 

the  new 

p » q • exp  (-y) 

(14) 

where 

9(0*^)  " ^1^10 

(15) 

Also  introduce  another  variable  r(y,*)  using 

s » r • exp  (-a2y) 

(16) 

where 

r(0,X)  ■ ^2^20 

(17) 

In  terms  of  the  new  set  of  dependent  variables  q,  r,  and  q and 
independent  variables  y and  \ the  rate  equations,  Eqs.  (11),  (12), 
and  (13)  reduce  to 

^ In  q • - aL  r]  (18) 

^ In  r • - a3  t)  (19) 
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- (1-n)  q e"bly  - Kj^  q q e”bly 

+ (l-q)r  e"b2y  - Kj  q r e'b2y  - q (20) 

where  the  constants 

\ - 1 + (a/Nc  oLl),  b2  - a2  + (P/N0  au)  (21) 

A FIRST  INTEGRATION 

An  important  integral  Is  obtained  by  eliminating  the  common  T)  - 
term  in  Eqs.  (18)  and  (19).  Integration  with  respect  to  the  y - 
variable  results  in 

r ■ C q i v 3 a^/a^  m ( c^j2  ’ ^1,2^ ^Ul  ” ^Ll^ 

As  a consequence  of  the  boundary  conditions  imposed  on  the  variables 
q and  r , the  arbitrary  function  of  integration  in  Eq.  (22)  is 
just  a constant  C.  Note  that  if  both  excitation  and  probe  beams  are 
of  different  wavelength  (implying  unequal  absorption  coefficients) 
one  expects  that  in  general  the  ratio  v ^ 1.  However,  for  equal 
wavelengths  the  ratio  is  identically  unity  (Ref.  11). 

The  three  first  order  partial  differential  equations  in 
Eqs.  (18)  through  (20)  are  reduced  to  the  pair  of  equations. 


In  q - - »i  n 

I?  * U'l)  q e’biy  ' K1  1 4 e’biy 

+ (1-q)  C qV  e b2y  - K2  q C qV  e b2y  - q 


The  two  level  system  is  now  described  in  terms  of  the  variables 
q,  r,  and  q,  governed  by  Eqs.  (22)  and  (23). 


REDUCTION  TO  A SECOND-ORDER  EQUATION 

For  purposes  of  numerical  evaluation,  it  is  more  convenient 
to  deal  with  a single  higher  order  differential  equation  rather 
than  with  a pair  of  first  order  equations.  Equations  (23)  are 
combined  to  form  a second  order  partial  differential  equation  in 
terms  of  one  of  the  variables,  q or  q.  By  far  the  simplest 
equation  is  obtained  if  the  q variable  is  eliminated  from  the 
pair  of  equations.  In  terms  of  a new  variable  z,  where 

In  q - z + b^  y (24) 

one  finally  has  the  lengthy  general  result 


This  second  order  partial  differential  equation  is  obviously 
"highly  nonlinear"  and  can  only  be  solved  numerically.  As  an 
alternative  we  will  consider  simplifications  that  are  practical 
as  well  as  physically  meaningful. 

Assume  that  interactions  with  the  probe  beam  can  be  neglected. 
That  is,  terms  in  Eq.  (4a)  containing  the  absorption  cross-sections 
aL2’  au2  are  to  be  neglected  in  this  analysis.  The  net  result 
is  to  discard  the  last  term  in  Eq.  (25).  (Surely  if  there  is  no 
probe  light  then  C » 0 and  this  term  will  vanish.)  With  this 
assumption  and  the  transformation 


» (y.*) 


In  (1+Kl)  + z 

ln  [ T <°L1  *L1  + °U1  V ] 


(26) 

+ In  Ix  (x,t) 
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then  Eq.  (25)  reduces  to 


- ' [(1+*“>  I?  + b3  e“  + bl] 
or,  dividing  by 


dO>jy)d>v 


<l«“>  S(b“ 


^ eu  + 1 

d(bjy)  bx 


where 


li  + (1+K1)  bl  . gUl^*Ll**Ul?  + — 


aLl*Ll+aUl*Ul  NoaLl 


(27) 


'3  ~ (1+KX) 

The  q-  variable  is  related  to  the  quantity  u as 

„ - (-bi/.l)  a + ^7) ) 

Equation  (27)  can  be  integrated  numerically,  subject  to  the 
initial-boundary  conditions 

u (0,70  - In  j\  (aTl*T1  + aTTl*,n)  | + In  A 


(28) 


(29) 


LI  LI  uUl  Ur 


] 


ilL 


^(bxy) 


* (ax  n (y,0)  + bL)/  bx  - - 1 (30) 


A-0 


where  a zero  initial  value  for  the  upper  level  population  has  been 
assumed  in  Eq.  (5)  (note  that  this  assumption  for  the  initial 
population  distribution  can  be  altered  without  modification  of  any 
prior  results).  Equation  (27)  is  recognized  as  a hyperbolic 
partial  differential  equation  in  canonical  form. 
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With  no  probe  beam  present  there  are  six  material  constant 
groups  that  must  be  known  in  order  to  achieve  a numerical  solution. 
Specifically,  these  quantities  are  b^,  K^,  and  a^  — terms  that 
appear  explicitly  in  Eqs.  (27)  and  (28).  In  addition,  the 
scaling  quantities  Nq  ol^  , t,  and  the  constant  term  in  the  boundary 
condition  of  Eq.  (30),  (namely,  *yi^  mU8t  *1®° 

given.  Furthermore,  if  a weak  probe  beam  is  present  in  the  system, 
then  a seventh  constant  must  be  known,  namely  v . Note  that  the 
quantities  and  A2  are  experimentally  controlled  parameters 
independent  of  the  material  understudy.  Likewise,  the  thickness 
of  the  sample,  x,  is  assumed  known.  In  the  following  sections,  we 
present  a few  experimental  procedures  that  can  be  used  to  determine 
some  of  these  needed  constants  from  a given  material  sample. 


I 
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EVALUATION  OF  THE  MATERIAL  CONSTANTS 

This  section  discusses  several  experimental  methods  that  are 
useful  for  material  parameter  estimation.  Specific  photochromic 
data  are  presented  as  a convenient  means  to  illustrate  the  procedures. 

STEADY  STATE  SOLUTION 

The  steady  state  conditons  are  found  from  Eq.  (27).  Setting 
the  partial  derivative  with  respect  to  equal  to  zero  and 
separating  variables  leads  to 


— I + eU  --  & _ - b 

(b1/b3)  + eu  dy  3 


(31) 


where  the  - symbol  denotes  the  steady  state  value  of  the  variable, 
(e.g.,  lim  u (y,*)  s u as  t -►••).  With  the  substitution 


w - (l>1/b3)  + eu 

one  has 


(32) 


L W"  (bi/b3) 

Upon  integration. 


u - (l-O*^))  • 


WbjL /b£" 
w 


dw 


bi  dy 


In 


u 


<bl/b3> 


] 


b^  y + const  (33) 


The  constant  of  integration  is  evaluated  using  the  boundary  conditon 
at  y ■ 0;  namely,  Eq.  (5c). 

The  variable  u (y,X)  and  the  originally  introduced  quantity  1^ 
(x,t)  are  related  by 
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u • In 


L T (<TL1  *L1  + aUl  *U1>  J + ln  rl  C34, 


Thus,  at  the  entrance  face  x — 0 the  steady  state  value 


i (0)  - In  Tt  (nL1  *L1  + a01  ♦U1)J  +ln 


+ In  A, 


Using  this  result  to  evaluate  the  constant  term  in  Eq.  (33),  one 
finally  has 


In  (I^/A^)  - I 1 - (bj/b3) 


- - b,  y - 1 - (b 


)j  In  |(I1/A1: 
<Vb3>)  lnV 


1 + b 


where  is  the  steady  state  value  of  1^  measured  after  the  excitation 
beam  has  traversed  a distance  x through  the  material.  The  constant 
b^  is  defined  as 

b5  = bl  J10/b3  A1  (1  + Kl>  (37) 

b.  (Vn  +J:} 

T “<c,Ll*Ll+0Ul*Ll)+No',Llc'ut(*Lr*Ul) 

The  quantities  b- , b3,  and  b$  are  known  in  terms  of  the  material 
constants  whereas  I and  x are  experimentally  controlled  or  measured 
values.  In  principle,  Eq.  (36)  can  be  used  to  help  evaluate  the  material 
constants  from  experimental  data.  Unfortunately,  the  functional 
form  of  this  equation  only  all  allows  for  "trial” and* error"  estimates. 


A useful  simplification  is  possible  if  we  consider  a "small 
signal  approximation"  where  (1*1)  for  which  i « 1.  The 

logarithmic  terms  in  Eq.  (36)  can  be  replaced  with  the  approximations 


In  (1-i)-  - i 

In  (l+b5-i/l+b5)  * - i/(l+b5) 
with  which  Eq.  (36)  simplifies  to 

b b (l+b5) 

i ^ — r — 7~r — — y » constant  * y (38a) 
bl  + b3  b5 

Thus,  for  the  "small  signal  case,"  the  fractional  change  in  transmitted 
excitation  light  is  linearly  proportional  to  the  thickness  y.  The 
slope  of  this  linear  relation  can  be  easily  related  to  several 
material  constants 


b,  b,  (l+bR) 


A,  a. 


bl  + b3  b5  " Ax  (1+Kl)  + J] 


+ 1 + 


No  °L1 


(38b) 


SIMULTANEOUS  TWO- BEAM  INTENSITY  MEASUREMENTS 

Figure  2 shows  the  experimental  system  used  to  measure  the 
transmitted  light  levels  in  both  excitation  and  probe  beams  after 
passage  through  a specific  photochromic  absorber  sample.  The  argon 
laser  448  nm  line  was  used  for  sample  excitation  and  the  helium-neon 
laser  633  nm  line  served  as  the  probe  light  in  this  demonstration. 

A photochromic  film  similar  to  that  employed  by  Baldwin  (Ref.  8) 

(type  63*071,  American  Cyanamid  Co.),  was  used  in  these  demonstrations. 
The  wavelength  sensitivity  of  this  material  is  well  suited  for  use 
with  the  blue-green  argon  laser  lines  for  excitation  together  with  the 
helium-neon  light  as  a probe  beam.  Both  laser  beams  were  spatially 
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filtered,  expanded,  and  collimated  to  illuminate  uniform1!'  * portion 
of  the  film.  A beam  splitter  directed  both  beams  onto  -ue  film 
and  shutters  were  used  to  control  exposure  times.  Transmitted 
beams  were  separated  using  narrow  band  pass  filters  and  detectors 
monitored  the  individual  changing  beam  intensities  in  time. 

At  every  instant  in  time  the  excitation  and  probe  beam 
intensities,  after  passing  through  the  absorber,  are  related  by 
Eq.  (22).  In  terms  of  the  measured  beam  intensities  one  derives  from 
Eq.  (22)  and  the  boundary  conditions  of  Eq.  (5)  the  result 

In  I2  - v In  1^  + [In  A2  * v In  A^  + y (v  b^  - b2)  ] (39) 

Thus,  a plot  of  the  quantities  In  I2  vs  In  1^  should  be  linear  with 
slope  v and  intercept  equal  to  the  value  of  the  last  three  bracketed 
terms  in  Eq.  (39).  Figure  3 illustrates  how  well  a straight  line 
represents  the  demonstration  experimental  results.  The  least 
square  fit  value  of  v » - 0.66  falls  within  Baldwin's  range  of 
values  (v  » - 0. 65  to  - 0.93).  No  numerical  estimates  were  obtained 
from  the  intercept  value  in  this  demonstration  since  the  experimental 
data  were  not  corrected  for  reflectivity  losses  at  the  front  face 
of  the  photochromic  film  and  only  relative  beam  intensities  were 
measured. 

POST- EXCITATION  DECAY 

The  time  constant  t can  be  found  from  decay  measurements  of  the 
excited  state  population  after  removal  of  the  excitation  beam.  The 
probe  beam  is  used  to  monitor  the  decay  rate.  We  first  assume  that 
the  interaction  between  probe  radiation  and  excited  state  population 
can  be  neglected  as  a consequence  of  using  a very  weak  probe  beam 
or  of  small  cross  section  values.  The  rate  equations  similar  to 
Eqs.  (19)  and  (20)  are  now 
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(40) 


_2_ 

dy 


In  r • • aj  T) 


(41) 


with  the  initial  distribution  t\  (y,0)=  (y)*0.  That  is,  an 

excited  state  population  is  created  and  exists  at  the  time  when 
the  excitation  beam  is  just  removed.  The  decrease  in  excited 
state  population  density  affects  the  probe  light  transmission 
character is  t ic  s . 

Integration  of  Eq.  (41)  is  straightforward,  resulting  in  the 
well-known  exponential  decay,  (y , X)  — q (y)  exp  (-X) . Sub- 

stituting this  result  into  Eq.  (40)  and  integrating  with  respect 
to  the  y-variable  finally  yields 

y 

nc  (y)  dy  (42) 


In  I2  » In  S2  + a3 


J 


The  quantity  S2  s I2  (y,0)  represents  the  initial  probe  light 
intensity  level  transmitted  through  the  absorber. 

The  change  of  the  helium-neon  probe  light  after  excitation 
of  the  photochromic  absorber  was  measured  in  the  demonstration 
experiment  previously  mentioned.  Figure  4 illustrates  the  variation 
of  the  probe  beam  light  level  at  several  points  in  time  after  the 
excitation  beam  is  removed.  In  accord  with  Eq.  (42)  the  quantity 
In  I2  should  vary  as  an  exponential  function  of  time.  The  curves 
shown  in  the  figure  represents  matches  to  the  experimental  data  pcints 
using  the  functional  form  of  Eq.  (42)  for  three  values  of  t.  A 
value  of  t - 1900  sec  was  found  to  offer  the  best  overall  curve 
fit.  This  estimated  value  for  t is  higher  than  that  reported  by 
Baldwin  (Ref.  8),  namely,  1100  < t < 1400  sec. 


A less  restrictive  model  allows  for  an  interaction  between  the 
probe  radiation  and  the  level  populations.  As  a practical  measure 
assume  now  that  the  interaction  involves  only  the  excited  level  by 
setting  0^2  ■ 0 in  Eqs.  (8)  and  (9).  Since  1^  (x,t)  - (x,t)  » 0 

after  excitation  one  has 


V ’ ' S3  J2 


-Py/n0  °L1 

T aU2  *U2  J2  116 


(43) 


The  quantity  J2  cannot  be  normalized  as  before  with  respect  to  the 
J2o”term  since  the  latter  is  not  defined  when  » 0,  (see  Eqs.  (10)). 
Instead,  we  normalize  with  respect  to  the  quantity  t au2  <|>u2 
by  defining 


S “ T °U2  *U2  J2 


(44) 


With  this  substitution  Eq.  (43)  becomes 


. o n C 

dy  a3  ^ S 


la  . - s n e 

d*  S 6 


-Py/N 


o Ll 


(45) 


- 


If  the  quantity  P is  zero  or  can  be  neglected  then  the  coupled 
Eqs.  (45)  can  be  integrated. 


With  the  transformations 


Q v 

S - e , tj  • e 


(46) 


Equations  (45)  finally  reduce  to 


M _ .a  g'v 
dy  a3  6 


& 
* d* 


- (1  + e^) 


(47) 


18 


Differentiation  of  the  first  of  Eqs.  (47)  with  respect  to  A 
yields 

(q  + eQ)  (48) 

Integration  first  with  respect  to  y and  then  with  respect  to  A 
results  in 


dQ 

Q+eQ+C 


- A + f (y) 


(49) 


where  C is  evaluated  from  the  boundary  conditions  at  y » 0 and 
is  given  by 


C » - In  [A2  t au2  4>u2  3 " A2  T au2  <I>U2  (50) 

Finally,  in  terms  of  the  variable  J2,  one  has 


(y,o) 


d J, 


J2  []  In  ( J2  M2  )+t  °u2*  U2  (J2_A2^!1 


-A+f(y) 

(51) 


The  integral  is  to  be  evaluated  numerically.  The  combination  of 
terms  t a^2  <t>u2  appear  as  a single  item  to  be  determined  utlimately 
from  the  numerical  integration  and  match  to  the  experimental  data. 
The  arbitrary  function  of  integration  f(y)  can  be  found  from  the 
steady-state  properties,  i.e.,  ^ (y,«>)  » 0 and  J2  (y,»)  — A2* 

Thus,  with  the  above  mentioned  approximations  the  material  constant 
t and  t au2  <*u2  can  be  found  from  experimental  data.  (Note  that 
Mohn' s (Ref.  9)  treatment  of  optical  bleaching  for  a photochromic 
material  suggests  an  approximation  for  which  Eq.  (51)  can  be 


integrated  to  estimate  the  material  constants.  When  applicable, 
this  estimated  value  could  be  used  in  a recursive  fashion  to  refine 
the  numerical  integration.) 

LINEARIZED  PARTIAL  DIFFERENTIAL  EQUATION  APPROXIMATION 

We  consider  simplifications  to  the  second  order  equation, 

Eq.  (27).  As  a first  assumption  set  a » 0.  Consider  an  expansion 
of  the  quantity  u about  a particular  quiescent  value  of  1^,  say 
< 1^  >,  in  the  form 

u » In  [t  (aL^  <J>L1  + <t>u]L)  < 1^  >]  + In  I^/<I^>  (52) 

using  Eq.  (26).  The  small  signal  approximation  Ij^  < 1^  > + 6 
suggests  that 

u » In  [t(^l1  ♦h  + aui  ♦yi^  <Il>  3 + ln  (1+6/<Ii>)  (53) 
~ uQ  + 6/<  lx  > 

where  u is  a constant.  With  this  simplification  Eq.  (27)  becomes 


Equation  (54)  is  a linear  partial  differential  equation  which  can 
be  readily  integrated  subject  to  the  initial-boundary  conditions 
of  Eq.  (30).  The  solution  for  this  linearized  problem,  expressed 
in  terms  of  the  quantity  1^  (x,t),  can  be  written  as 
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In  1^  InA,  - 


h 

biL 


l+(l-W1)<I1>/J10 


+ 1 


(55) 


ai<1 £/jip 


+ l+(l+K1)<I1>/J10  y exp 


•( 


(1+K-i  )<I-i  > 

1 + * I X 

J10 


(Note,  the  second  term  in  Eq.  (55)  would  be  identical  to  a "small 
signal"  steady  state  result  of  Eq.  (38)  in  accord  with  the  assumptions 
<I^>  ^ J^q  ~ A^.  The  last  term  in  Eq.  (55)  vanishes  in  steady  state.) 

Figure  5 shows  several  measured  1^  - t data  points  from  the 
demonstration  experiment  and  a curve  fit  determined  from  the 
functional  form  of  Eq.  (55)  as  A + B exp  (-Ct).  The  constants  A, 

B,  and  C were  evaluated  from  the  measurements  at  the  later  time 
intervals  t » 400,  500,  and  600  sec.  From  the  computed  exponential 
factor,  C,  and  using  the  value  t •=“  1900  sec  obtained  previously, 
one  find  that  the  quantity  (1+K^  <Ij>/j 3.  For  example, 
assuming  the  value  ~ 1 and  <I^>  of  the  order  of  1 p,  watt/cm^, 
one  finds  that  J^q  is  also  of  the  order  of  1 p.  watt/cm^.  (Baldwin's 
(Ref.  8)  equivalent  value  for  J^q  is  about  20  p watts/cm^.)  In 
addition,  the  computed  value  of  B in  the  above  mentioned  exponential 
form  leads  to  an  estimate  of  the  product  a^  y ^ - 10.  This  result 
Implies  that  the  ratio  aui/aLl  muc^  less  than  unity  as 
expected  (see  Eq.  (10)). 
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The  accomplishments  of  this  work  are  threefold:  First,  general 
rate  equations  have  been  developed  either  as  a pair  of  first  order 
or  a single  second  order  partial  differential  equations.  As 
previously  noted,  these  equations  have  a broad  scope  and  are 
useful  to  model  absorption  phenomena  and  energy  level  population 
changes.  The  procedure  was  shown  to  be  valuable  for  the 
identification  of  significant  material  constant  groupings  as 
characteristic  constants;  that  is,  several  scale  parameters  and  non- 
dimensional  groups  are  new  physical  insights  to  the  problem. 

Secondly,  the  reduction  of  the  descriptive  equations  to  a hyperbolic 
partial  differential  equation  in  canonical  form  is  another  Important 
result  that  implies  relative  simplicity  and  accuracy  for  numerical 
evaluations.  In  addition,  an  integral  solution  was  derived 
relating  probe  and  exictation  beam  intensities  in  the  form  of  a time- 
independent  power  law  (see  Eqs.  (22)  or  (39)).  This  result  suggests 
a new  experimental  technique  for  monitoring  the  excitation  beam. 
Finally,  several  experimental  methods  were  presented  in  which  the 
calibration  test  data  could  be  analyzed  to  estimate  a number  of 
material  constants.  A photochromic  was  used  as  a convenient 
demonstration  media. 
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Optical  System  to  Measure  Attenuation  of  Collimated 
Excitation  (488  nm)  and  Probe  (633  nm)  Beams  Traversing 
a Two- Level  Photochromic  Absorber 


TIME  t (SEC) 


Fig.  4 Measured  Changes  of  the  Probe  Light  Level 
Transmitted  Through  the  Absorber  after 
Excitation  and  Several  Comparison  Exponential 
Time  Functions. 
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TIME  t (SEC) 


Fig.  5 A Comparison  of  Measured  Excitation  Light  Levels 
Transmitted  Through  the  Absorber  with  a Linear 
Approximation  Solution 


APPENDIX  A - NUMERICAL  INTEGRATION  PROGRAM 

To  accommodate  a numerical  integration  procedure,  the 
hyperbolic  differential  equation  of  Eq.  (27)  can  be  expressed  as 


" £<u'  ilu/a<V»*  f 


(A-l) 


together  with  the  initial  value-boundary  value  conditions  of  Eq.  (30) 
written  as 


u (o,X)  - C 


(A-2) 


_ du 
d(bjy) 


at  A ■ 0 or  equivalently 


(y,o)  - - bxy  +C 


(A- 3) 


For  convenience  set  b^  • 1,  (equivalent  to  redefining  the  y-varlable). 

Integration  of  Eq.  (A-l)  in  terms  of  the  differential  or  step 
size  A*  yields 


- & 


+ f AX 


(A-4) 


• y , X-»A  * 


Also  the  change  in  the  u-variable  for  a step  change  y can  be 
written  as 

u[(y+Ay),  (*HA?oQ-  u(y,(7v+A?0)  + -Ay  (A-5) 

y,*+A> 


I 

l 


To  obtain  the  recurrence  formulas  upon  which  the  integration 
scheme  is  based  let 

X - n • AX  , y - s • Ay  (A-6) 

where  n,  s are  Integers.  Thus,  Eqs.  (A-4)  and  A-5)  become 


Is,  n+1 

u ((s+1),  (n+1))  - 


+ f(s,n)  AX 

s,n 

u (s,(ntl))  + |j} 


(A- 7) 

Ay  (A-8) 


1 s,n+l 

together  with  the  initial  value-boundary  value  conditions 
u (o,n)  - C 


• - 1 or  u(s  ,o)  » - 8 • A y + C 


8,0 


(A-9) 


For  simplicity  in  notation  define 


°.,n  ‘ “<»•"> 


D.  - 

s,n  dy 


s,n 


Finally,  the  recurrence  formulas  are 


(A-10) 


Ds,n+lV*Ds,n  + f <8'n*X 
Us+l,n+l  " Us,n+1  + Ds,n+1  * A y 


(A-ll) 


Uft  - C 
o,n 


u - - s * Ay  + C 

8,0 

A-2 


A program  for  the  solution  of  Eqs.  (A-ll)  is  given  in  Figure 

A-l.  The  Interpretations  of  several  program  variable  are  stated 

in  Table  A-l.  Initial- boundary  conditions  are  stored  on  files, 

(see  program  statement  1242).  The  values  of  U and  D 

8 p n 8 p Q 

computed  during  the  run  are  also  stored  on  files,  (see  program 
statement  1320).  Solution  of  the  recurrence  equations  is  carried 
out  between  statements  1372  and  2220.  An  output  file,  (see  statement 
2050)  consists  of  exponentiated  u-values  at  a fixed  value  of  n(orA) 
over  a range  of  s(or  y)  values. 

At  the  starting  value  of  the  interger  n,  the  change  in  the 
variable  U8  with  s is  calculated  over  the  range  of  s-values  and 
subsequently  stored.  The  value  of  u is  increased  by  unity  and  the 


computations  are  repeated.  The  final  computed  values  of  U and 

s p n 

D are  stored  in  file  format,  to  be  used  as  input  values  when  the 

o | n 

program  is  run  again. 


Note  that  in  this  version  of  the  program  only  two  steps  in  n are 
utilized  due  to  computer  storage  limitations,  i.e. , N 9*2.  To  go 
beyond  two  steps  requires  program  changes  in  the  arguments  of 
statement  numbers  2140  and  3600. 


TABLE  A-l  - SIGNIFICANT  PROGRAM  VARIABLES 


Computer 

Program 

Variable 


C 

B 

T 

L 

M2 

N9 

M5 

SI 

Nl,K 

F 

U(sl,K) 

D(slfK) 

S 

G3,G4 


Mathematical 

Equivalent 

lnCT(aLl*Ll  + -oiWJ  +lnAl 

Vbl 

A* 

Ay 

smax 

n 

max 

Initial  value  of  n 
s 
n 
f 

Us,n 

Ds,n 

exp  (U) / exp  (C) 
accuracy  test  variables 


A- 4 


PHOTOS 

1000  LONG  U7?  D7»  S 

1010  LONG  DC  1500*  2 ]j  U[  1500*  2 3 

1020  LONG  C » B j R ? B4  j T > L > M5  > M2  ;•  G 1 » G2  j F 

1021  LONG  FI 

1030  DIM  HtC  8 3 

1031  DIM  G*[8  3>F«8  3 
1040  DIM  U$[  8 ]> DtC  8 3 
1050  DIM  E$[8  3»H$C8  3 
1060  MRT  D=ZER 

1070  MRT  U=ZER 
1080  FILES  *)*)*?*>* 

1090  PRINT  "INPUT  BOUNDARY  CONDITION  C=?" 

1100  INPUT  C 

1118  PRINT  "INPUT  PRRRMETER  B«?" 

1120  INPUT  B 

1130  PRINT  "TIME  STEP  T=?  AND  SPRCE  STEP  L=?" 

1140  INPUT  T 

1141  INPUT  L 

1150  PRINT  "INPUT  THE  NUMBER  OF  STEPS  IN  L M2=" 

1160  INPUT  M2 

1170  PRINT  "INPUT  THE  STRRTING  TIME  STEP  M5=" 

1180  INPUT  M5 

1185  PRINT 

1186  PRINT  "HOW  MANY  TIME  STEPS  N9?" 

1187  INPUT  N9 
1190  C1=C 
1200  61=B 
1210  T1=T 
1220  L1=L 
1230  M3=M5 
1240  M4=M2 

1242  PRINT  "LIHAT  ARE  THE  INITIAL  CONDITION  FILE  NAME 

1243  INPUT  Ft 

1244  INPUT  Gt 

1245  PRINT  Ft  n Gt 

1246  ASSIGN  Ft»3»Yl 
1248  ASSIGN  Gt>4»Yl 

1250  PRINT  "BOUNDARY  CONDITION  C="5C1 
1260  PRINT  " PARAMETER  IN  DIFF.  EQ.  B="5B1 
1270  PRINT  "TIME  STEP  T="5T1 
1280  PRINT  "SPACE  STEP  L="5L1 
1290  PRINT  "STARTING  TIME  STEP  M5="?M3 
1300  PRINT  "NO.  OF  SPACE  STEPS  M2="5M4 


Fig.  A-l  Computer  Program  for  the  Solution 

of  the  Recurrence  Formulas  of  Eq.  (A* 11) 
(1  of  3) 


1301  PRINT 

1302  PRINT  " H9=“,H9 

1310  PRINT  “START  NEXT  TIME  STEP  JJ5  = "5M5+N9 

1311  PRINT 

1320  PRINT  " INPUT  U*  AND  b*  FILE  NAMES" 

1321  PRINT 

1 o30  INPUT  U* 

1331  INPUT  D*  ' 

1340  ASSIGN  U$ > 1 , Y 1 
1350  PRINT  Y1S 
1 360  FibS  I IjN  D$  ’ 2 , V 1 

1370  PRINT  Y1»U$»D$ 

1371  PRINT 

1372  FOR  R7= 1 TO  N9 

1380  09=0 

1390  Sl  = l 
1395  K=0 

1400  FOR  Nl*l  TO  M5+R7+1 

1410  READ  #3! FI 

1412  IF  K=2  THEN  1430 

1415  IF  H 1 < ( M5+R7-2 ) THEN  1430 

1417  K=K+ 1 

1420  UC  SI , K ]=F 1 

1430  NEXT  HI 

1435  RESTORE  #3 

1440  FOR  Sl-1  TO  012-1 > 

1450  READ  #25 D7 

1460  DC  SI  1 3=H7 

1470  NEXT  SI 
1480  FOR  S 1 = 1 TO  M2 
1490  READ  #15U7 

1500  UE  SI , 1 ]=U7 

1510  NEXT  SI 
1520  K= 1 

1530  FOR  Sl=l  TO  012-1) 

1540  F= 1 +EXP < UE S 1 j K ] ) 

1550  F = < F * D C S 1 , K ] ) + < E:  * E X P < U E S 1 , K ] > ) + 1 

1560  F=-F 

1 570  DE  S 1 1*  K+ 1 1— DC  S 1 , K ]+F+T 

1 588  UE S 1 + 1, K+ 1 3=UE S 1 j K+ 1 ]+HE  S 1 , K+ 1 3*L 

1590  G1=DE  SI , K+l  T~ DE 31 j K ! 

1600  G3=.0001 

1610  IF  G 1 >G3  THEN  1650 

1 620  G2=UC  S 1 + 1,  K + l 3— UE  S 1 . r.+ 1 3 

1630  G 4 =.@001 

1 640  IF  G2\=»G4  THEN  1690 

1650  PRINT  "ABOVE  THRESHOLD" 

1660  09=1 

1670  PRINT  C, B, T, L, M5, M2, SI 

1680  81012 

169 0 NEXT  31 

1700  IF  09=1  THEN  2220 

1850  K=2 


Fig.  A-l  Computer  Program  for  the  Solution 

of  the  Recurrence  Formulas  of  Eq.  (A-ll) 
(2  of  3) 

A- 6 


I860 

RESTORE  #1  - 

.13-70 

RESTORE  452  \ .. 

4900 

FOR  J=1  TO  M2 

1910 

A=UC J»K3 

1920 

PRINT  #15  A 

1930 

NEXT  J 

1935 

GOTO  3000 

1940 

FOR  J=1  TO  CM 

1950 

B4=DC  J»  K 3 

I960 

PRINT  #25 B4 

1970 

NEXT  J 

1980 

PRINT  #15  END 

1 990 

PRINT  #2 5 END 

2000 

RESTORE  #1 

2010 

RESTORE  #2 

£040  NEXT  R7 

2050  PRINT  "INPUT  THE  OUTPUT  FILE  HR ME  H$= 
£060  INPUT  H* 

2061  PRINT 

£070  CREfiTE  Y1 » H$  > 2 » £56 

£080  IF  Y 1=0  THEN  £110 

£090  PRINT  " CAN’T  OPEN  THE  H*  FILE" 

£100  GOTO  £220 
2110  ASSIGN  H$ > 5 » Y 1 
£120  PRINT  Y15H* 

2130  FOR  J-l  TO  M2  STEP  <M£/38> 

2140  S=UC  J j 2 ]-C 
2150  S=EXPCS> 

£160  PRINT  #55 S 

2170  PRINT  S 

£180  NEXT  J 
2190  PRINT  #5 5 END 
2210  RESTORE  #5 
2215  PRINT  #4  5 END 
2220  END 
3000  S1=M2 

3130  FOR  J=1  TO  CM5+R7-2) 

3200  READ  #45  FI 

3500  NEXT  J 

3600  F1=UC  SI > 2 I 

3700  PRINT  #45  FI 

3705  IF  R7=N9  THEN  4000 

3800  RESTORE  #4 

4000  GOTO  1940 
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